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ABSTRACT 

We perform a detailed study of the bispectrum of the Sunyaev-Zel'dovich effect. Using an analytical 
model for the pressure profiles of the intracluster medium, we demonstrate the SZ bispectrum to be 
a sensitive probe of the amplitude of the matter power spectrum parameter as- We find that the 
bispectrum amplitude scales as -Btsz oc crl l ~ 12 , compared to that of the power spectrum, which scales 



as A t sz oc er 8 ~ . We show that the SZ bispectrum is principally sourccd by massive clusters at 
redshifts around z ~ 0.4, which have been well-studied observationally. This is in contrast to the SZ 
power spectrum, which receives a significant contribution from less- well understood low-mass and high- 
redshift groups and clusters. Therefore, the amplitude of the bispectrum at £ ~ 3000 is less sensitive to 
astrophysical uncertainties than the SZ power spectrum. We show that current high resolution CMB 
experiments should be able to detect the SZ bispectrum amplitude with high significance, in part due 
to the low contamination from extra-galactic foregrounds. A combination of the SZ bispectrum and 
the power spectrum can sharpen the measurements of thermal and kinetic SZ components and help 
distinguish cosmological and astrophysical information from high-resolution CMB maps. 
Subject headings: cosmology: dark matter — galaxies: clusters: general — intergalactic medium 



1. INTRODUCTION 

The Sunyaev-Zel'doy ich (SZ) effect 

(|Sunvaev fc Zeldovichl 119801 ) is the principal source 
of CMB temperature anisotropy power at angular scales 
smaller than a few arcminutes. It comprises of two 
components, known as the 'thermal' and 'kinetic' effects. 
The former arises from inverse Compton scattering of 
CMB photons off hot electrons pervading galaxy cluster 
environments, distorting the Planckian form of the 
CMB spectrum. The kinetic SZ effect is caused by 
the Doppler shifting of CMB photons via Thomson 
scattering off clouds of electrons with a non-zero bulk 
velocity (along the line-of-sight) relative to the CMB rest 
frame. Recently, both the South Pole Telescope (SPT 
Lueker et all 12010: IShirokoff et 10120111 IReichardt et~ail 



theoretical uncertainty on the amplitude of the tSZ 
signal (for a fixed cosmological model). Precise pre- 
dictions of the SZ power spectrum have been ham- 
pered by the uncertainties in modeling the state of 
the intra-cluster medium (ICM) over a broad range 
of mass and redshift , and at large cl uster-centric radii 
(jBattaglia et al .1120101: [Trac et al.ll20TTl S10). Early mod- 
els and simulations produced predictions for the ampli- 
tude of the thermal SZ power spectrum that are dis- 
crepant with recent measur ements by more than a fac- 
tor of 2 (|Lueker et al]|2010D. Recent work, using hy dro 
dynamical simulations (jBattaglia et al.l I201O | 2011[). N 



body simulations+semi-analvti c gas models flTrac et al 



2011 



2011 ) and analytic models (SlO.IEfstathiou fc Migliaccio 



: iChaudhuri fe Mai umdar 2011), have significantly 



20111) and the Atacama Cosmology Telescope (ACT, 



Fowler et af1l20lQlDunkfev et al.H201lD have placed con- 



straints on the amplitudes of the kinetic and thermal SZ 
contributions t o the small-scale CMB power spectrum. 
For instance, IReichardt et al.l (J2011I ) measured the 
thermal SZ amplitude at 150 Ghz to be 3.69 ± 0.65^K 2 
and placed an upper limit on the kSZ amplitude of 
< 2.8/xK 2 (at 95% confidence and assuming no SZ-point 
source correlation). 

Measurements of the thermal SZ power spectrum 
probe the abundance-weighted pressure profiles of hot 
gas in groups and clusters over a wide range of red- 
shifts. This is interesting in a cosmological context 
because the dependence on the abundance of clusters 

produces a sensitiv e scali ng with as\ i.e., Cj oc g« 
(iKomatsu fc Seliakl [20021: iTrac et all [20111 iShaw et al.l 
120101 denoted as S10 hereafter). However, efforts to 
constrain as in this way have been inhibited by a large 



reduced the tension between the observed and predicted 
values. However the distribution of amplitudes between 
different models and simulations is s till significantly 
large r than the measurement errors ([Reichardt et al] 
120 111 ), degrading the constraints that can be placed on 
cosmological parameters. Furthermore, as current exper- 
iments do not cover a sufficiently broad frequency range 
to enable the tSZ and kSZ signals to be distinguished 
spectrally, the kSZ contribution must currently be mod- 
eled and subtracted from the total SZ power, providing 
an additional source of systematic uncertainty. 

In this work we propose the SZ bispectrum as a new, 
robust method for constraining the thermal SZ ampli- 
tude. The bispectrum, or its real-space counterpart, the 
3-point function, is widely used to measure the devia- 
tion of a signal from Gaussianity both in the CMB and 
in large scale structure data. The bispectrum in the 
CMB map can be measured by using three CMB tern- 
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perature values and taking the average over the survey 
sky. The primordial CMB bispcctrum can be used to 
search for signatu res of primordial non -Gaussian fluc- 
tuations (see e.g.. iKomatsu et al.1 120111 and references 
therein). The dominant secondary anisotropy signals in 
the CMB temperature power spectrum arise from lensing 
of the CMB by large-scale structure, the ISW effect and 
the SZ effect. These signals, either alone or from their 
cross correlation, give rise to a measura b le bispectrum 
(ICooravl 120001 : IKomatsu fc Soergell 127)01 ICoorav fc Hul 
l2000h . The full non-Gaussian information (primordial or 
from low-redshift structure) can be obtained by Fourier 
transforming the CMB temperature field and taking the 
product of three temperature values that form a triangle, 
then summing over all possible triangle configurations in 
harmonic space, assuming all harmonic modes forming 
the triangles arc independent. 

One can also measure a skewness spectrum in which 
two sides of the triangle are summed over to mea- 
sure the skewness as a function of the third side 
(|Munshi fc Heavens! I2010T ). The total skewness spec- 
trum in a CMB map consists of that generated by CMB 
lensing, the thermal and the kinetic SZ effects, point 
sources, a signal arising due to the cross correlation of 
these components, and any primordial contribution. A 
significant fraction of the kinetic SZ power spectrum is 
generat ed by Gaussian d ensity fluctuations in the linear 
regime (jShaw et al|[2011l ). Although patchy reionization 
may c ontribute a significan t fraction of the kSZ spec- 
trum (Mesing cr et al.M2012h . it is reasonable to expect 
that the kinetic SZ bispectrum is very small - the overall 
kSZ amplitude is 50-75% smaller than the tSZ. The SZ 
skewness spectrum is thus likely to be dominated by the 
thermal SZ effect. As the primary CMB power decreases 
rapidly at £ > 2000, the CMB lensi ng bispectrum should 
also becomes negligible at high £. ICoorav et al.l (J2000I ) 
have shown that the SZ skewness can be measured us- 
ing the Planck data with a S/N«10, assuming a perfect 
subtraction of the primary CMB and other foregrounds. 

In this work, we focus on the SZ bispectrum, showing 
that it can provide information complementary to that 
from the SZ power spectrum. First, as the SZ skew- 
ness spectrum is dominated by the thermal SZ compo- 
nent it does not suffer from the observed thermal-kinetic 
SZ degeneracy. A combination of the SZ skewness and 
power spectra could therefore help disentangle the ther- 
mal and the kinetic signals. We show that the data from 
current small-scale CMB experiments have the potential 
to constrain the thermal SZ amplitude using bispectrum 
measurements with 5 — 10% accuracy. Furthermore, a 
tightened constraint on the kinetic SZ amplitude can also 
improve our understanding of the reionization scenario 
(IShaw et alJ[20TTI : IZahn et alll20ll . 

Second, as we will show, the signal is dominated by 
massive clusters at intermediate redshift for which high- 
precision X-ray observations exist. This is in contrast 
to the power spectrum where the signal mainly comes 
from the l ower mass and hi gher redshift groups and clus- 
ters (e.g.. lTrac et al.ll2011l ). The theoretical uncertainty 
in the SZ skewness spectrum is thus expected to be sig- 
nificantly smaller than that of the SZ power spectrum. 
Combined measurements of the power spectrum and the 
bispectrum can thus be used to distinguish the contribu- 
tion to the power spectrum from different cluster mass 



and redshift ranges. 

The goal of this study is to derive the SZ skewness 
spectrum using an analytic model for the ICM gas and 
analyze its astrophysics and cosmology dependence. In 
Section 2 wc discuss the theory of the SZ bispectrum, 
and in Section 3 wc describe our model for calculating 
the thermal pressure of the intra-cluster medium used 
in the bispectrum calculation. In Section 4 we present 
our results, including the distribution of the skewness 
spectrum signal over mass and redshift, the theoretical 
uncertainty, cosmological dependence, and the prospects 
of detectability from the current data from SPT like 
surveys. Section 5 presents summary and discussions. 
Throughout this work, we assume a fiducial cosmology 
with zero curvature and parameters: Qb = 0.045, S7 m = 
0.27 0-8 = 0.8, h = 0.71, n s = 0.97, (jKomatsu et all 
l20Tl and w = -1. 



2. THEORETICAL FRAMEWORK 
2.1. Sunyaev-Zel'dovich effect 
The Compton y-parameter for the thermal SZ effect is 
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(1) 



where I is the distance along the line-of-sight; n e and 
T e are the electron number density and temperature, re- 
spectively. 

The thermal SZ power spectrum can be calculated by 
summing up the square of Fourier transform of the pro- 
jected SZ profile, weighted by the number density of clus- 
ters of a given mass and redshift, 



Ci = f{x„f 



, d V f ,, „rdn(M,z)_. 2 .,._ „, 



.(2) 
where V(z) is the comoving volume per steradian, 
f(x u ) = x v (coth(x I ,/2) - 4) with x v = hv / '(/cbTcmb) 
at frequency v and is negative (positive) below (above) 
220 GHz, Tcmb is the temperature of CMB at z = an d 
n{M, z) is the halo mass function of lTinker et al.l (|2008f) . 
y(M,z,£) is the Fourier transform of the projected SZ 
profile of a cluster, given by 



y(M,z,£) 



(3) 
where x = r/r s , £ s = DA{z)/r s , r s is the scale radius of 
the 3D pressure profile, Da(z) is the angular diameter 
distance at redshift z and P e is the electron pressure pro- 
file. Note that we use the mass function fit for 400 times 
the mean density, whereas the pressure profiles are de - 
fined at the virial overdensity (jBrvan fc Normanlll998lh 
We convert the mass function definition to the virial den- 
sity in Eq. [2] The cosmology dependence of the SZ power 
spectrum comes from the angular diameter distance and 
the growth function while the astrophysical dependence 
comes from the pressure profiles. 

2.2. Bispectrum Theory 

The CMB temperature fluctuations A8(h) in a certain 
direction, n, can be expanded into spherical harmonics 
as 
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The angular bispectrum is 
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where the angle-averaged quantity in the full sky limit 
can be written as 
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which has to satisfy the conditions: m\ + ra2 + ma = 0, 
4 + 4 + 4= even, and % - £ j \ < 4 < ^ + 4- Here 
5(444) is the bispectrum in the full-sky limit. In the 
small angle limit, the flat-sky approximation is valid. By 
defining 6(444) to be the bispectrum in the flat-sky 
limit, then the correspond ence betw een the full-sky and 
the flat-sky bispectrum is (Hu 2000), 



£(444) < 



(24 + 1) (24 + 1) (24 + 1) 
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Throughout this work, we focus on the high-^ regime 
(to avoid confusion from primary CMB and its lensing), 
where a flat sky approximation is valid. We will refer 
to the flat-sky bispectrum; namely, 6(444) as the bis- 
pectrum for brevity and .5(444) as the full-sky bispec- 
trum. Note that, for t > 500, the Wigner 3-j symbol can 
be approximated as 



4 4 4 




(_1)V2 



7r[Z(L-24)(£-24)(£-24)]V4 

(8) 
if L = 4 + £2 + £3 is even, and vanishes for odd L. The 
approximate expression for Eq. [5] is valid to better than 
1 percent for £ > 500, the £ range we are interested in 
this study. 

The thermal SZ bispectrum is the volume integral of 
the cube of the Fourier transform of the pressure profile 
weighted by the halo mass function, 



6(444) = f{x v ) 
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x y(M, z, ti)y(M, z, h)y(M, z, 4)- (9) 

We sum over the two smaller values of £s and writing the 
bispectrum as a function of the largest £, we define the 
skewness spectrum as 



A(4 = ./E 62 («^), 



(10) 



where £ + 4 + 4 is even and 4 < 4 < (■■ 

In previous work, the amplitude of the thermal SZ 
power spectrum has been defined in terms of the ex- 
pected amplitude at I = 3000 for a given set of cos- 
molo gical parameters ([Lueker et al.ll2010t Shiroko ff et al.l 
[2WllReichardt et aljfcoilUDunkley et al.ll2011l ). In this 
study, we define A^sz as the ratio of the thermal SZ power 
spectrum amplitude at non-fiducial cosmology to that at 
og = 0.8 as 

. _ £3000(0-8) ni s 

tsz = r TZ nsV <• > 

£3000 (,o- 8 = 0.8) 



Parameters 


Fiducial 


Range 


Dark Matter concentration (Ac) 


1 


0.8- 1.2 


Energy Feedback (ej) 


4 x icr 7 


(l-lOO)xlO- 7 


Dynamical Friction Heating (erjjw) 


0.05 


0.0-0.1 


Non-thermal pressure normalization (00) 


0.2 


0.0-0.3 


Non-thermal pressure evolution (/3) 


0.5 


-1-+1 



TABLE 1 
Fiducial values and the uncertainty range in the ICM 

PARAMETERS. 

Similarly, we define the amplitude of the skewness spec- 
trum as 

n _ A 300 o(o'8) / 10 n 

-tftsz = -7 -, -r-zz, (U) 

A3000(CT8 = U.Bj 

where A3000 is the amplitude of the skewness spectrum at 
£ = 3000. For the fiducial cosmology, B t sz = -Atsz = 1- 
Given that the skewness spectrum is proportional to the 
cube of the pressure profile the expected relation between 
B^z and A t sz is -Btsz c< Alg Z , assuming similar mass 
and redshift distribution of the power spectrum and the 
bispectrum. In Section 14.31 we study the relation be- 
tween the SZ amplitude of the power spectrum (A t sz) 
and the skewness spectrum (-Btsz) in more detail. 

3. MODELING THE INTRA-CLUSTER MEDIUM 

To calculate the thermal pressure profiles of groups 
and clusters, we use the physically motivated analytic 
model presented in S10. Here we briefly describe the key 
features of the model, and refer the readers to S10 for 
more details. 

The model assumes that intra-cluster gas resides in hy- 
drostatic equilibrium within the gravitational potential 
of the host dark matter halo, 



dPtotjr) 
dr 



-Pg( r ) 



d$(r) 
dr 



(13) 



where p g (r) is the gas density at radius r from the clus- 
ter center, the total pressure, Ptot(^) = Pth(r) + -fnt(r), 
is given by a sum of thermal and non-thermal pressures, 
and <j)(r) is the gravitational potential. The gas is as- 
sumed to have a poly tropic equation of state, Ptot = 
Po(p g /po) r , where T = 1.2 and Pq and po are the central 
gas pressure and density respectively. The total mass dis- 
tribution is ass umed to follow the Navarr o-Frenk- White 
(NFW) profile (|Navarro et al.lfl996l [l997h . 



Ptot(r) 



x(l + xf 



(14) 



where p s is the normalization constant, x = r/r s , and 
r s is the characteristic radius and can be defined in 
terms of the concentration of halos, r s = r v j r /c v j r where 
"vir" refers to the virial overdensity (|Brvan fc Normanl 
I199&T ) with respect to the critical density of the uni- 
ve rse. We adopt the concentration-mass relation given 
in lDuffvet aPpoMsh : 



c vir (M, z) = 7.85^ c 



M m 



2 x lQ w h- x M G 



-0.081 



(l+z)-°- 7 \ 
(15) 
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Fig. 1, — Comparison between the observed 'universal' pressure profile (scaled with respect to P500) of (Arnaud ct al. 20f3) (blue dashed 
line, shaded region denotes 20% scatter) and our model predictions at z=0 for a cluster of hydrostatic mass Ml^ E = 3 X 10 14 /i -1 Mq. In 
each panel we vary one of the input model parameters to demonstrate their effect on the pressure profile. The black solid line represent 
the fiducial values as given in Table [T] The top-left panel shows the profiles obtained when ef is 10~ 5 (red dotted line) and ejjM varies 
from 0.2 to 0.0 (green dashed lines, from top to bottom at 2.R500); in the top right panel we vary A c from 0.8, 0.9, 1.1, 1.2 (red dotted 
lines, top to bottom at 2.R500), the bottom-left panel shows the variation with cxq from 0.3, 0.24, 0.14, 0.0 (red dotted lines, bottom to 
top at 2.R500) and the bottom-right panel shows the variation with /3 at z — 0.5 over {1, 0.5, 0, —1} where self-similar evolution is factored 
out (red dotted lines, from bottom to top). Note that the observed result is at < -R50O1 the blue dotted line is an extension of the Arnaud 
profile based on the simulation result. 



where A c is the normalization factor, with the fiducial 
value A c = 1. Note that the gas is distributed as a 
massless tracer of dark matter distribution. 

The model accounts for non-gravitational processes 
that affect the thermal properties of the ICM including 
star formation, energy feedback from supernovae (SNc) 
and Active Galactic Nuclei (AGN), and non-thermal 
pressure support due to bulk gas motions and turbulence. 
Star-formation is implemented by converting a certain 
(mass and redshift-dependent) fraction of the gas mass 
into stars. We adopt the observed stellar mass fraction 
for t he local X-ray groups and clusters (jGiodini et all 
l2009f ) and assume that the stell ar populations ev o lve ac - 
cording to the "fossil" model of iNagamine et aLl (|2006l ). 



Following Os triker et all (|2005l ). a fraction of the rest 
mass energy of the stars, e/M*c 2 , is put back into the 
gas due to energy feedback from SNe and AGN, where 
ef is a free parameter. The model also accounts for 
the transfer of energy from dark matter to gas by dy- 
namical heating by infalling substructures, determined 
by the product of the total binding energy of the halo 
I-^dmI and a free parameter £dm- Finally, we take into 
account the non-thermal pressure support due to ran- 
dom gas motions as seen in hydrodynamical simulations 
(e.g-lKav et alll200l iRasia et alJ[200l lLau et ai]|2059t 
iBattaglia et alj2010l;lNelson et al.lboill ). We assume the 



non-thermal pressure to be a certain fraction of the total 
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Fig. 2. — Contribution to the SZ skewness spectrum of objects 
within a given mass range (top) and redshift range (bottom) range. 
We show the ratio of the skewness spectrum in a mass/redshift 
range to the total skewness spectrum. 



pressure and a power law in radius, 
-ftot V -K500 



(16) 



where a(z) = aof(z), olq is the ratio of non-thermal to 
total pressure at R500, the radius at which the spher- 
ical overdensity of the cluster is 500 times the critical 
density of the universe (enclosing mass Afsoo). f{z) = 
min[(l + z) 13 , (4~™ nt / a o - 1) tanh(/3z) + 1], n nt and f3 are 
the radial and redshift dependence of the non-thermal 
pressure support, respectively. We set n nt = 0.8, mo- 
tivat ed by comparisons to hydrodynamical simulations 
(SlO.rBattagli aet al.ll2011[ ). leaving j3 as a free parame- 
ter. 

In summary, our ICM model has five free parameters 
Ac, e/, £dm, «0j and /3. Since the key observable for the 
SZ effect is the pressure profile, we use mea surements 
of the pressure profiles bv lArnaud et al.l (|2010T l based on 
31 massi ve (M 5n o > lO^h^M ^), low redshift (z < 0.2) 
clusters. lArnaud et al.l (|2010[) snowed that the pressure 
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Fig. 3. — Top panel: the distribution of the skewness spectrum 
over mass and redshift at £ = 3000. At M 500 = 5 X 1O 13 /i -1 M , 
the lines show the cumulative distribution function, CDF(>M, 
<z)=0.1 to 0.9 from bottom to top with a step of 0.1. Bottom 
panel: the distribution of the SZ power spectrum at £ = 3000 (for 
comparison with the skewness distribution). The lines represent 
similar quantities as that in the top panel. 



profiles in their sample adheres to a universal form with 
^2 0% scatter. Fi g ure [1 ] shows the pressure profile fit to 
the lArnaud et a l. (2010) data along with our predictions 
for a typical cluster of mass M^ S E = 3 x 10 14 h- 1 M Q 
at z=0, where M| u q E is a cluster mass estimated assum- 
ing the hydrostatic equilibrium (see S10 for more discus- 
sions). Each of the 5 parameters are varied over a wide 
range, including some extreme cases with a very high 
feedback parameter or zero non-thermal pressure. 

The parameter e/ determines the amount of non- 
gravitational ('feedback') energy injected into the ICM. 
High levels of feedback inflates the gas distribution, low- 
ering the pressure in the inner region while boosting 
it in the cluster outskirts. For cluster-mass halos, the 
feedback energy is a small fraction of the binding en- 
ergy and has a minimal effect on the pressure profile. 
The upper- left panel of Figure [1] shows how e/ changes 
the pressure profile in a cluster-mass halo at z = 0. 
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Increasing e/ by factor of ~ 100 changes the pressure 
profile by few percent at radii R < -R500 and increases 
to about 30% at 3i?soo- This is in qualitative agree- 
ment with the findings from hydrodynamical simulations 
including AGN feedback which indicate the impact of 
AGN fe edback is more pronounced in groups than in 
clusters flBhattacharva et al.ll2008tlPuchwein et al.ll2008l; 
Siiacki et al.ll2007t iTevssier et al.H201lHMcCarthv et all 



2010HBattaglia et al.ll2010Ll201lHShort et al.ll2012lT ~We 
note however that our feedback model is highly simplis- 
tic and does not capture the full detail of AGN feedback 
as implemented in hydrodynamical simulations (for ex- 
ample, energy injection via jets). 

£dm determines the fraction of the total dark matter 
energy transferred to the ICM during mergers. Clusters 
frequently undergo mergers and hence €dm may have 
a substantial effect on the ICM physics of clusters. As 
shown in Figure Q] (top- left), €dm changes the pressure 
profile by 20-30% over the entire radius range as £dm 
varies from — 0.2. 

The normalization of the concentration-mass relation 
(A c ) determines the characteristic radius of an NFW 
dark matter profile, approximately corresponding to the 
radius at which the slope changes from —1 to —3. The 
upper-right panel of Figure [T] shows the change in the 
thermal pressure profile as A c varies from 0.8 — 1.2. In- 
creasing the concentration deepens the central potential 
which draws the gas towards the cluster core, steepening 
the pressure profile. The impact of increasing halo con- 
centration provides a similar (although opposite) effect 
to that of energy feedback. 

The nonthermal pressure parameter, ao, determines 
the fraction of the total gas pressure contributed by ran- 
dom gas motions. As ao is varied is the thermal pres- 
sure fraction adjusts accordingly such that the total gas 
pressure is fixed. As shown in Figure [1] (bottom- left), 
the thermal pressure profile changes by about 60% at 
i?50o as ao varies from — 0.3. /3 determines the redshift 
evolution of the non-thermal pressure fraction. S10 ob- 
tained an estimate for f3 by co mparing the redsh ift and 
1 outputs of the simulations of lLau et al.1 (J2009I ) , finding 
that the best fit value for ao increases from 0.18 — 0.26 
implying f3 = 0.5. The bottom-right panel in Figure Q] 
shows the pressure profile of a cluster at z=0.5 and how f3 
changes the thermal pressure profile. A positive value of 
(3 increases the amount of non-thermal pressure at higher 
redshift and hence decreases the thermal pressure as red- 
shift increases. The pressure profile decreases by about 
20% as /3 varies from —1 to 1. 

Overall, we find that the range e/ < 10~ 6 , eoM = 
0.04 - 0.06, A c = 0.8 - 1.1, and a = 0.14 - 0.24 
to be in agreement with the Arnaud pressure profile 
within the uncertainty (see Figure Q}. Note that, as 
described in S10, the constraints on e/ < 10~ 6 and 
£dai = 0.04 — 0. 06 mainly come from t he gas fraction- 
mass relation of IVikhlinin et al.1 (J2006I ) . We also find 
the entropy scaling relation (e.g., entropy - X-ray tem- 
perature relation) predicte d by our ICM mod el to be in 
agreement with the data of lPratt et al.l (|2010f) within the 
ranges of the ICM parameters. This illustrates that the 
current observations can indeed constrain the ICM pa- 
rameter space significantly. 

4. THE SZ SKEWNESS SPECTRUM 



Mass Range 
MsooIH^Mq] 


power 

% 


z-range 


power 

% 


< 10 14 


3 


0-0.3 


38 


(1 -2.5) x 10 14 


15 


0.3-0.6 


45 


(2.5 - 5) x 10 14 


30 


0.6-0.9 


12 


(5-25) x 10 14 


52 


> 0.9 


5 



TABLE 2 

Distribution of the skewness spectrum over different 

mass range at £ ~ 3000. 

4.1. Mass and redshift contributions 

Figure [5] shows the fraction of skewness spectrum sig- 
nal contributed by different ranges of halo mass and red- 
shift. At £ ~ 3000, galaxy groups (M 500 < 1O 14 /i _1 M ) 
source only 5% of the signal, while lower mass clus- 
ters contribute about 20%. The majority of the sig- 
nal comes from the massive clusters. The distribution 
of the skewness spectrum is tabulated in Table [5] over 
different mass and redshift ranges. We also calculate 
the signal expected from clusters that are below the de- 
tection threshol d of the curr e nt ge n eration of SZ sur- 
veys. Following IBenson et al.l (|2011[ ) ; IVanderlinde et all 
(|2010t) . we choose the SPT mass threshold of M500 ~ 
2.5 x 10 14 /i _1 Mq, roughly corresponding to the detec- 
tion significance of 5cr. The skewness spectrum gets only 
20% (at £ = 1500) to 40% (at £ = 6000) of the signal from 
clusters with M500 < 2.5 x 10 H~ 1 Mq. In contrast, the 
SZ power spectrum gets about two-thirds of the signal 
from groups and clusters with M500 < 2.5 x 1O 14 /i -1 M0. 

Figure [31 top panel shows the cumulative fractional 
distribution of the skewness spectrum for halos with mass 
< M500 and redshift > z, from 10-90% of the signal (top- 
left to bottom-right) at £ = 3000. In terms of the redshift 
distribution, about 90% of the skewness spectrum signal 
comes from z< 1. The bottom panel of Figure [3] shows 
the distribution of the SZ power spectrum, showing that 
a significant fraction of the thermal SZ power comes from 
high-z, low-mass objects (i.e., about 40% of the power 
coming from z > 0.7 and M500 < 1O 14 /i _1 M ) where the 
ICM properties arc still very uncertain. In comparison, 
the SZ skewness spectrum signal is dominated by low-z, 
massive clusters. 

4.2. Astrophysical Uncertainties 

We now evaluate how robust the predicted skewness 
spectrum is to uncertainties in the underlying ICM model 
parameters. Figured] shows the variation of the SZ skew- 
ness spectrum over the range of ICM parameters given 
in Table [1] and plotted in Figure [TJ As discussed in Sec- 
tionO €f does not have a significant effect on the massive 
clusters that dominate the skewness spectrum. Conse- 
quently, as shown in the top-left panel of Figure [TJ the 
skewness spectrum varies only by 14% as e/ varies from 
(1-100) xl0~ 7 . On the other hand, the pressure profiles 
of massive clusters vary substantially with £dm , and the 
skewness spectrum varies by about 40% as ei)M is varied 
from to 0.2. In contrast, the power spectrum receives 
a significant contribution from a wider range in mass, 
M 500 = 10 13 - 10 15 /i _1 M o , especially at high £. As a re- 
sult, e/ shows more variations at £ — 8000 (group scales) 
than at £ = 2000 (cluster scales, see Fig. 7 of S10). 
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Fig. 4. — Dependence of the SZ skewness spectrum on the cluster gas physics. The solid line indicates the fiducial model with ef = 
4 X 10 — 7 ,£dm = 0.05, A c = 1.0,oo = 0.2, and /3 = 0.5. The dotted and dashed lines represent the skewness spectrum predicted by our 
ICM model for a range of ICM parameter values. The top-left panel shows the impact of tf (red dotted line: 10~ 5 and black solid line: 
4 X 10 — 7 ) and eoM (green dashed lines: 0.0, 0.04, 0.06, and 0.2 from top to bottom); the top-right panel shows the change with A c 
(0.8, 0.9, 1.0, 1.1, 1.2 from bottom to top); the bottom-left panel shows the variations with «o (0.0, 0.14, 0.2, 0.24, 0.3 from top to bottom); 
and the bottom-right panel shows the variation with /3 (—1,0,0.5,1 from top to bottom). 



Increasing the normalization of the concentration-mass 
relation (A c ) deepens the potential well of the clusters 
which steepens the pressure profiles. The top-right panel 
of Figure [4] shows the change in the skewness spectrum 
is about ± 30% for a ± 20% change in A c (compared to 
the power spectrum where the variation is ~ 20%). This 
can be understood from the fact that the skewness ampli- 
tude is oc(power spectrum) 15 , so a 20% change in power 
spectrum corresponds to a 30% change in the skewness 
spectrum. Also note that in contrast to the case of ey, 
the effect of A c is not significantly mass-dependent. 

The lower-left panel shows the effect on the skewness 
spectrum of varying the non-thermal pressure parameter 
«o- The skewness spectrum changes by about 70% across 
the range o.q = — 0.3. As the SZ skewness spectrum gets 
most of its signal from massive clusters at z < 0.4, the 



uncertainty in the z-evolution of the non-thermal pres- 
sure (/3) contributes only about 20% uncertainty to the 
SZ skewness spectrum. The SZ power spectrum, on the 
other hand gets its signal from high-z objects and as a 
result (3 adds ~ 30 — 40% uncertainty in the SZ power 
spectrum (bottom-right panel of Figure 7 in S10). 

Figure [5] shows the astrophysical uncertainties in the 
shaded area around the fiducial model. Here, the total 
uncertainty is computed by adding the uncertainties (cal- 
ibrated by the observed pressure profile and discussed in 
detail in Section [3]) in each ICM parameter in quadra- 
ture. On average, about 75 — 80% of the skewness spec- 
trum signal arises from massive lower redshift objects; 
i.e., M 500 > 5 x l0 14 h~ 1 M Q at (z)0.4, and hence is less 
sensitive to the variation in ICM parameters than the 
power spectrum (c.f. Figure 7 in S10). Figure [5] shows 
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Fig. 5. — Astrophysical uncertainties in the SZ skewness spec- 
trum. The solid line represents the fiducial model and the shaded 
region indicates astrophysical uncertainty for the ICM parameter 
range calibrated by the observed pressure profile (see Section [3j . 
At I = 3000, the astrophysical uncertainties are ~ 33%. 
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the gas physics uncertainty in the skewness spectrum to 
be about 33%. We also provide the theoretical predic- 
tion of the skewness spectrum for the fiducial values of 
cosmological and ICM parameters. 

At this point, it is worth commenting on effects of sev- 
eral simplifying assumptions in our model. First, we as- 
sumed a constant V = 1.2. However, recent work sug- 
gests that r could be radially dependent, especially at 
large radii (r < i? V ir)- However, the bulk («95%) of the 
skewness signal arises from r < R y\ r , and T ~ 1.14 — 1. 3 
over the redshift range z = — 1 (|Battaglia et al.ll201lD . 
We therefore find that the skewness spectrum changes 
by about ±10 — 15% in the specified range of I\ Also, 
there is uncertainty in the number of baryons converted 
into stars in clusters, since this is only constrained by 
redshift zero observations. Secondly, the redshift depen- 
dence of the model assumes that evolution of the stellar 
fraction is the same for all galaxies, although it is pos- 
sible that the star fraction evolve differently in a cluster 
environment. In order to assess the impact of star for- 
mation model on the skewness spectrum, we considered 
a no evolution model and found that the difference be- 
tween the "no-evolution" mo del and the fiducial mo del 
(namely the "fossil" model of iNagamine et aTl (J2006T )) is 
negligible. Finally, recent work on the thermal SZ power 
spectrum have shown that variations about a mean pro- 
file may produ ce up to a 15% excess in the spectrum at 
£ ~ 3000 (e.g., iBattaglia et alJl201lh . This variation is 
non-Gaussian and radially-dependent and can lead to the 
enhancement of the bispectrum amplitude. 

4.3. The relation between the power spectrum and the 
skewness spectrum amplitude 

Given a measurement of the SZ skewness spectrum am- 
plitude, wc can predict the expected amplitude of the 
thermal SZ power spectrum, namely the -Btsz — -Atsz 
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£ 2 /(27r)A(0yuK 3 


1000 


0.66 


2000 


0.69 


3000 


0.56 


4000 


0.44 


5000 


0.35 


6000 


0.28 



TABLE 3 

The theoretical prediction of the skewness spectrum for 

the fiducial values of icm and cosmology parameters. 



relation. Recall that the power spectrum and the skew- 
ness spectrum are proportional to the square and the 
cube of the pressure profile, respectively. If we assume 
the skewness spectrum and the power spectrum have the 
same mass and redshift distribution, then any changes 
in gas physics affects the skewness spectrum and (power 
spectrum) 15 in the same way. However, as shown in 
Sec 14.11 the skewness spectrum depends mostly on the 
massive clusters at intermediate redshift, while the power 
spectrum receives a substantial contribution from higher 
redshift galaxy groups. Furthermore, the cluster physics 
affect the gas properties in lower mass structures more 
than in clusters. Therefore, the -Btsz — -<4tsz relation is 
expected to deviate from the above relation and depends 
on both cosmology and cluster astrophysics. 

To investigate the dependence of the Btsz — ^tsz rela- 
tion on cluster physics and cosmology, we vary erg from 
0.7 — 0.9 in steps of 0.05 and compute the -Btsz — ^tsz 
relation for each value of a% for our fiducial gas physics 
model. Figure [S] shows that the skewness spectrum am- 
plitude is proportional to the power spectrum amplitude, 



Btsz(^ = 300 °) = Km = 3000) 



(17) 



This is shallower than the expected ^4tsz behavior, pri- 
marily because the skewness spectrum signal comes from 
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more massive clusters which are not point-like at the an- 
gular scales equivalent to t ~ 3000. We thus study the ra- 
dial dependence of the skewness spectrum signal and find 
that less than 5% signal comes from [l-2]i? v ; r , 20% comes 
from [0.5-l]i? vir , 60% from [0.2-0.5]i? vir , the rest 15% 
from < 0.2i? v ; r . Thus a simple scaling of the pressure 
normalization with mass cannot explain the i?tsz — -Atsz 
relation. Next we investigate the variation of the scal- 
ing relation with the ICM parameters. Within the range 
of the ICM parameters allowed by the observed pressure 
profile, we find the Btsz — -<4tsz relation is robust and 
the overall uncertainty does not vary by more than 7% 
(indicated by shaded area in Figure [5]). The tight scaling 
arises because both the power spectrum and the skewness 
spectrum depend on the pressure profile. As a result, the 
amplitudes of the two spectra varies in similar ways as 
we vary our model parameters. 

To check the robustness of the relation further, we con- 
sider two extreme gas physics scenarios, such that the 
SZ power spectrum varies from 2.5 — 10/^K 2 (roughly 
by factor of ~ 4) which is roughly twice the theoreti- 
cal uncertainty assumed in the curre nt SZ power spec- 
trum analysis (jReichardt et al.ll201lD . The first case is 



where the pressure profile amplitude is maximally sup- 
pressed by high feedback and non-thermal pressure and 
low concentration: e/ = 10~ 5 , euu = 0.1, ao = 0.3, 
and A c = 0.8. The second case other case consists of 
the parameter choice that predicts a higher amplitude: 
e/ = eon = «o = 0.0 and A c = 1.2. We find that even 
for the two extreme gas physics scenarios the Asz — -Btsz 
has an overall uncertainty of only by ~ 15% (indicated by 
dashed lines in Fig[6]). The slope of the relation changes 
from 1.33 to 1.47 between these two cases. This is en- 
couraging as using the skewness spectrum data we can 
measure i?tsz , then use this theoretical relation to derive 

Atsz- 
The relation in Eq. [17] implies the skewness spectrum 

amplitude varies as er^ 1-12 . We also investigate how the 
skewness spectrum varies with other cosmological pa- 
rameters. Figure [7] shows the variation of the skewness 
spectrum with a% and fib- Each panel shows how each 
cosmological parameter changes the skewness spectrum 
while other parameters are fixed at their respective fidu- 
cial values. The figure shows that the skewness spectrum 
depends very sensitively on as, and the amplitude varies 
as Bg Z oc ct| 16 , as expected. The variations with Hi, is 
moderate, while ft m , h, wq, and n s show much smaller 
variation. We find the SZ amplitude for the skewness 
spectrum follows the scaling (at t = 3000), 
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(18) 



0.97 



ft, 



0.27 



-0.46 



Note that not all cosmological parameters follow a simple 
-Btsz = ^tsz relation. The difference in scaling is due to 
the different range of masses in the mass function con- 
tributing to the skewness spectrum and the power spec- 
trum. Note that we ignore the cos mology dependence 
of the concentration-mass relation (jBhattacharva et al.1 
l2011af) in our study which can lead to a slightly different 
scaling of the skewness spectrum amplitude with cosmol- 
ogy. We provide the theoretical prediction of the skew- 
ness spectrum in Table 3 for our fiducial cosmology and 
gas parameters. 

5. DETECTABILITY 

We now assess the detectability of the SZ skew- 
ness spectrum in the presence of other secondary CMB 
anisotropy signals and extra-galactic foregrounds. Other 
than the thermal SZ effect, point sources (dusty, star- 
forming galaxies (DSFGs), and radio sources) will also 
contribute to the skewness spectrum signal. The main 
difference is while the thermal SZ effect contribute a neg- 
ative b ispectrum, the point sources cont ribute a positive 
signal (Rubino-Martin & Sunyacv 2003) . Assuming that 
the clustered component of the DSFGs follows the Gaus- 
sian distribution, the skewness spectrum of the corre- 
lated DSFG is zero. As mentioned before, the kinetic SZ 
follow approximately Gaussian distribution and is thus 
assumed to have negligible skewness. We also ignore a 
possible skewness spectrum signal arising from the cor- 
relation between the SZ effect and the DSFG sources. 

5.1. Noise Estimates 
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Fig. 8. — Detectability of the SZ skewness spectrum at 150 GHz. The left panel shows the SZ (black solid) and the point sources- 
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signal- to- noise integrated to a certain £ of the SZ skewness spectrum (black solid), the total- SZ+DSFG+radio (red dotted). We assume 
an SPT-like survey with sky coverage of 2500 deg , 18 /^K-arcmin, 1.2 arcmin at 150 GHz. 



The noise of the total skewness spectrum, assuming the 
primary contribution to the variance of the bispectrum 
co mes from the power spectru m of the CMB sky, is given 
bv lKomatsu fc Spergell (pOOlh 



iV 2 (£i£ 2 £ 3 ) = C(h)C(l 2 )C(£ 3 )A M2i3 



(19) 



where C{£) is the total power spectrum that includes con- 
tributions from (1) the SZ power spectrum computed us- 
ing Eq.[2] (2) the lensed CMB power spectru m computed 
using CAMB for the fiducial cosmology (jLewis et al.1 
l2000f K (3) the noise due to the finite resolution of an 
experiment, and (4) point sources, including dusty star- 
forming galaxies (DSFG) and radio sources. Note that 
Af 1 £ 2 £ 3 = 6 if all three £s are equal, Af^ 2 ^ 3 =2 if two 
of the three £s are equal and Ag ± g 2 g 3 =1 if all three £s 
are different. We describe the estimate of (3) and (4) in 
more detail below. 

The nois e due to the finite resolution of the ex periment 
is given by iKnoxl (|1995[ ); Uungman et al.l (|1996l ) 



N b (£) 



r £202 



w 



cxp 



8 log 2 



(20) 



where w~ 1 = [a viyi 9 /Tqmb] 2 is the weight per unit solid 
angle, <7 p i x is the noise per pixel, 9 is the FWHM of 
the instrument in radians, Tcmb is the CMB tempera- 
ture. Eg. 1191 gives the expression for noise in the full sky 
limit. For the fractional coverage, Eq. [19] is multiplied 

3/2 

by / skv , where / s k y is the fraction of the sky covered 
by the assumed surveys. Here we consider an SPT like 
survey with / s k y given by 2500 deg sky coverage, cr p i x 
is 18 /iK-arcmin at 150 GHz with the pixel size set by 
a Gaussian beam of FWHM=1.2 arcmin. The powers of 
the unresolved DSFG sources at 150 GHz are given by 
5.25 x 10" VK 2 (Poisson) and 4.4 x 10- 6 (£/3000)- 1 - 2 /zK 2 
(clustered), and that of the radio sources is 9 x 10~ 7 /iK 2 



(Poisson) at 150 GHz (|Reichardt et alj|20ll . 



5.2. Point Source Bispectrum 

The DSFG population comprises of a clustered compo- 
nent and a Poisson component, both of which contribute 
to the point source power spectrum. We assume that 
the clustered component produces Gaussian fluctuations 
and hence does not contribute to the bispectrum of the 
point sources. We thus consider the bispectrum contri- 
bution only from the Poisson component of the DSFGs 
and the radio galaxies. The Poisson component of the 
DSFG typically comprises of two populations: (1) the 
bright poin t sources which are detected by the CMB ex- 
periments ([Vieira et al.ll2010t IMarriage et al.|[20TH) and 
(2) the faint sources which follow a shallower distribu- 
tion compared to the bright sources, and are undetected 
by the current CMB experiments but contribute a signif- 
icant portion of the arcminute scale DSFG power mea- 
sured in the CMB experiments. Both populations are 
expected to contribute to the total point source bispec- 
trum. 

Given the DSFG flux number counts dN/dS, we can 
write the bispectrum of DSFGs as 
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(21) 

where the flux number counts is assumed to vary 
as dN/dS ex S~^ -1 , 5 m i n is the flux cut above 
which all sources are detected and removed from the 
CMB maps, -/Vdsfg(> Smin) is the number density 
(per unit solid angle) of the DSFG with S > S m - m , 
g {x v ) = 1/(65.55 MJysr" 1 )[sinh(x„/2)/a; 2 ] 2 , and x v = 
v/(56. 8GHz). Similarly, we define the bispectrum of the 
radio sources in terms of the radio source power spectrum 
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We assume all point sources above 5a at 150 GHz are 
detected and removed (masked). For the S PT 150 GHz 
band , this corresponds to S' m i n =6.4 mJy (|Vieira et al.l 
120101) . We adopt the number density of t he DSFG to be 

Nbsfg = 0.3 deg" 2 (jNegrello et al.ll2007t ) and /3 D sfg = 
1.93 for the total (bright +faint) DSFG population to 
be consistent with the total DSF G power spectrum mea - 
sured in the CMB surveys (e.g.. iReichardt et al.ll20il"l ). 
For the radio population we adopt A rac i(> S mnl ) = 
1.29 dcg~ 2 and /3 rac i = 1.03, which is consisten t with 
the 150 GHz source counts in iVieira et al.1 (J2010I ). The 
skewness spectra of the radio sources and DSFGs are cal- 
culated the same way as the SZ skewness spectrum by 
plugging Eqs. OH] and OH in Eq. [TO] 

The skewness spectrum for the various components at 
150 GHz are shown in the left panel of Figurc[8] Provided 
the point sources that are detected above 5a in the SPT 
band-power are removed, the DSFG and radio popula- 
tion make negligible contribution to the skewness spec- 
trum at £ < 4000. Beyond this scale the skewness spec- 
tra of point sources and SZ effect become comparable. 
If no point sources are removed, then the skewness spec- 
trum of the point sources (DSFG+radio) is a factor of 4-5 
higher compared to the case with point source removal 
and becomes the dominant signal at £ > 5000. Since 
the point sources and the SZ bispectrum have opposite 
signs, the total skewness spectrum (SZ+point sources) is 
smaller than the SZ skewness spectrum alone, especially 
at higher-^ where the point source bispectrum becomes 
non- negligible. 

5.3. Signal-to-Noise 

Finally, using Eq.0 we dehne the signal-to-noise (S/N) 
in terms of th e full-sky bispectrum. The S/N integ rated 
to a certain £ (jKomatsu fc Spergelll2001l: lHvill2000l ). A(< 
£) or Xi can be written as 



\(<e) 
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v B 2 (£ 1 £ 2 £ 3 ) 
„4l ^ 2 (W 3 ) : 



(23) 



where Is < I2 < h < I- The right panel of Figure |5] 
shows the total S/N for the SZ skewness spectrum. For 
a 2500 deg 2 SPT like survey at 150 GHz, the integrated 
S/N for the SZ skewness spectrum is ~ 16. Note that 
the total (SZ+point source) skewness spectrum is lower 
than the SZ skewness spectrum alone, because these two 
spectra have opposite sign. We find the S/N of the total 
skewness spectrum is only 10% smaller than the SZ only 
case. 

6. SUMMARY AND DISCUSSIONS 

We investigated the SZ bispectrum and the prospects 
of its detection with the current and upcoming CMB 
surveys. We define a skewness spectrum of the SZ effect 
which is a sum over all possible triangle configurations 
over the two smaller sides and expressed as a function of 
the largest side. About half of the SZ skewness spectrum 
signal arises from massive (M 500 > 6.5 x 1O 14 /i _1 M ) 



clusters in the local universe (z < 0.4). This is in con- 
trast to the SZ power spectrum whose signals come pri- 
marily from high-z, low-mass galaxy groups (M500 < 
2 x 1O u /i _1 M0 and z > 0.6) that have little observa- 
tional constraints at present. The skewness spectrum, 
therefore, is less susceptible to the large astrophysical 
uncertainties associated with high-z, low-mass groups. 

Our model for the SZ profile from groups and clus- 
ters reproduces the pressure profile of massive X-ray 
clusters in the local universe. This provides significant 
constraints on the ICM models of the present-day mas- 
sive clusters leaving only the redshift dependent param- 
eters unknown. Adopting wide priors on these param- 
eters, we determine astrophysical uncertainties in the 
skewness spectrum to be about 33%. This is in con- 
trast to the SZ power spectrum which is uncertain at 
the level of 50%, primarily because a significant frac- 
tion of the power comes from high-z, low- mass groups 
with Af 500 < 2 x lO^h-HlQ and z > 0.6. Such ob- 
jects have not been studied in as much detail as their 
high luass counterparts, making their theoretical model- 
ing and hence the SZ power spectrum highly uncertain. 
Other than the gas physics uncertainties, the theoretical 
modeling of the cluster mass function also add to the un- 
certainties in the bispectrum prediction. Because of the 
exponential nature of the mass function at the massive 
cluster end, a small change in cluster masses cause sig- 
nificant difference in the a bundance of massive clusters 
(jBhattacharva et al.|[20ll . 

The relation of the amplitude of the SZ skewness spec- 
trum and the SZ power spectrum is robust to these as- 
trophysical uncertainties. Indeed, the overall uncertainty 
of the relation is only 7% within the allowed range of 
ICM parameters. The skewness and the power spec- 
trum amplitudes of the thermal SZ effect is related as 
-Bg Z = ^4tszi consequently, the skewness spectrum varies 
as Bg Z oc erg 1 ' 6 . Given the astrophysical uncertainty of 
33%, the skewness spectrum measurements have the po- 
tential to constrain og to better than 3% accuracy. Since 
the SZ skewness spectrum gets its contribution primarily 
from the thermal component, the SZ skewness spectrum 
measurements can constrain the thermal SZ component 
uniquely We can then use the skewness spectrum mea- 
surements to break the degeneracy between thermal and 
the kinetic power spectrum measurements. 

We show that the current high resolution experiments, 
such as ACT, Planck, and SPT, have the potential to 
measure the skewness spectrum. For example, using the 
full 2500 deg 2 of data, SPT should be able to detect the 
SZ skewness spectrum with a signal-to-noise ~ 16. In- 
deed, the SZ skewness signal we predict in this study is 
detec ted by the ACT collaboration at 5a (|Wilson et al.l 
I2012T) . We note however that our noise calculation does 
not account for all possible sources, as we make several 
simplified assumptions. For example, the contribution of 
the trispectrum is assumed to be negligible and we ignore 
any spatial correlation between the tSZ and CIB. How- 
ever, as most of the SZ trispectrum signal comes from 
massive clusters (jShaw et al.l 120091) . masking out these 
objects should help reduce the SZ trispectrum contribu- 
tion to the non-Gaussian component of the variance of 
the skewness spectrum. Also we assume the lensing 3- 
pt function and its correlation with the SZ 3-pt function 
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is zero, because CMB lcnsing dominates at I < 2000 
while the SZ skewness spectrum gets most of the signal 
at £ ~ 4000 (see Fig 8, right panel). In this study we also 
assume the kSZ bispectrum is negligible. While this as- 
sumption is valid for the linear part of the kSZ effect, the 
non-linear kSZ effect contributes a non-zero bispectrum 
which is not included in our calculation. However, given 
that the overall magnitude of the kSZ effect is about 20- 
40% of the thermal SZ amplitude, the non-linear kSZ 
bispectrum should be small. Nevertheless, in order to 
use the SZ skewness spectrum as a robust cosmological 
probe, these assumptions must be tested using mock cat- 
alogs. We will present such study in our future work. 
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